High-accuracy transmission and fluorescence XAFS of zinc at 10 K, 50 K, 100 K and 150 K using the hybrid technique

High-accuracy measurement of the X-ray absorption fine structure (XAFS) of zinc metal in a temperature series from 10 K to 150 K, including 298 K, simultaneously in transmission and fluorescence using the hybrid technique is presented: the first transition metal XAFS using the hybrid technique and at low temperatures; and the first (hybrid-like) experiment at the Australian Synchrotron. A methodology for relative measurements and a methodology for cryostat measurements, and an approach to normalization, calibration of transmission measurements for scattering and other systematics to reference room-temperature data and yield a defined uncertainty are presented.


Introduction
X-ray absorption spectroscopy (XAS) has been a major field in physics for more than 100 years, and thousands of XAS papers are published each year. It enables investigation of widely used physical parameters such as the mass attenuation coefficient, atomic form factors and nanostructure in both solids and solutions. The development of the X-ray extended range technique (XERT) by Chantler et al. (2012a) has enabled a greatly improved accuracy to be achieved in X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS), with improvements often by orders of magnitude, and through this improvement discrepancies with theoretical predictions have been revealed. Despite the high dependence on XAS techniques in a wide range of scientific fields, very few experiments using standard X-ray absorption fine structure (XAFS) configurations have sufficient accuracy to test the method and current theoretical understanding. The development of the hybrid technique (Chantler et al., 2012bIslam et al., 2015Islam et al., , 2016Schalken & Chantler, 2018;Trevorah et al., 2019) is an attempt to overome the inadequacies of conventional XAFS systems. The hybrid technique allows for simultaneous measurement of transmission and fluorescence spectra. These are usually measured independently, but here can be directly compared. It has been found valuable for: dilute solutions including ferrocene and derivatives; for determining conformations when all shells of the different conformers are at the same radii; for nickel complex solutions; and for distinguishing local geometries where coordination numbers are identical. It has not, however, so far been applied to (solid or transition) metal foils. Also, it has not previously been applied to temperature series or low-temperature measurements.
Low-temperature measurements, usually with only a single small sample in the cryostat, also usually involve samples which cannot have their total surface profile mapped in the beam. A detailed two-dimensional xy scan is important to measure the local thickness in the beam and the integrated column density [t] and hence for absolute measurement of the attenuation coefficient ½=. This lends itself to the idea of using the hybrid technique and the use of some external reference (a secondary reference foil outside the cryostat which can be mapped), where normalization can proceed via transmission measurement. Here we use the same samples for fluorescence and transmission (i.e. internal and external references) to calibrate several possible systematics, in order to attain higher accuracy and insight.
Here we investigate the thermal evolution of zinc metal foil. Zinc is fairly well studied by XANES (Kraft et al., 1996), though repeated studies have found energy discrepancies of 2 eV or so (Eisa et al., 2005). A temperature series of X-ray diffraction (XRD) (Nuss et al., 2010) has provided valuable data on changes of the separations of the lattice sites. Several studies have investigated XAFS around zinc sites in compounds or complex systems (Gilbert et al., 2002;Chung et al., 2000;Sier et al., 2020). An XAFS study of zinc metal (Rae et al., 2010) yielded too sparse a grid for detailed analysis, though a recent work at room temperature using XERT (Ekanayake et al., 2021a,b;Sier et al., 2022) laid a strong and detailed baseline from which thermal evolution could be investigated. This is important for reference nanostructure and cross-platform portability when comparing analysis from XERT and hybrid, fundamental studies of zinc, and also for studies of small-scale atomic displacements and thermal expansion, the robustness of XAFS software, and particularly the extraction and understanding of the inelastic mean free path of the photoelectron and associated plasmon structure (Bourke & Chantler, 2014).

Hybrid attenuation measurements: experimental
The experiment was performed at the Australian Synchrotron XAS beamline using the hybrid technique (Chantler et al., 2012bBest & Chantler, 2022) to obtain XAFS data for zinc metal foils at temperatures of 10 K, 50 K, 100 K and 150 K. Zinc foils with nominal thickness of 10 mm and 50 mm were sourced from Goodfellow with reported purities of 99.9% and 99.95%, respectively. The foils were rolled, which can provide some static disorder from twinning or stacking faults, compared with an ideal crystal structure for example, but the foils exhibited no obvious pattern or defects. Mounts had low stress so had no impact on the metal nanostructure. They were placed inside a cryostat at 45 from the incident beam with a 100-pixel monolithic germanium detector, the fluorescence detector, placed at 90 (Fig. 1). Three ion chambers were utilized in the setup with one upstream (I 0 ), one downstream (I 1 ) and a third used to collect reference foil information with reference foils situated between I 1 and I 2 for energy calibration. A daisy wheel with different aperture sizes and range of aluminium foils was placed between the sample and each of the first two ion chambers. Fig. 2 presents the ratio of upstream (I 0 ) and downstream (I 1 ) ion chamber readings. Each point of the spectrum was measured three times providing a point-wise statistical uncertainty for each ion chamber reading. The high quality and resolution of the spectra are immediately visible by the detailed oscillations, low uncertainty and height of the first peak relative to the background. Fig. 3 reports the counting rates and statistics for the individual ion chambers. Foil thicknesses are chosen according to the expanded Nordfor's criterion (Nordfors, 1960;Chantler et al., 2000aChantler et al., ,b, 2012a with electronic settings optimized for counting statistics. At the highest level of absorption, the 10 mm foil maintains highly detailed structure and statistics, whilst the 50 mm foil goes below the noise floor. The upstream ion chamber shows a very smooth response function, with some Bragg peaks (Bragg glitches) that are characteristic of the Si(111) double-crystal monochromator used to monochromate the desired photon energy. Due to the high linearity of the ion chambers, these Bragg peaks disappear when normalized with the upstream counts. Fig. 2 shows the result of taking the ion chamber ratio for each temperature. Increased broadening for higher temperatures is already visible. Raw counts and experimental standard errors of three repeated measurements of the upstream and downstream ion chambers. Temperatures are given for the (a) 10 mm and (b) 50 mm foils. The 50 mm foil attenuated excessively towards the dark-current limit. Pre-edge counts read 1.3 Â 10 5 and 3.5 Â 10 4 for the 10 mm and 50 mm foils, respectively. (c) I 0 is independent of temperature and remained relatively consistent with the exception of Bragg glitches.

Figure 2
Transmission ratio of the downstream (I 1 ) ion chamber normalized by the upstream (I 0 ) ion chamber for (a, b) the 10 mm foil and (c) the 50 mm foil. The 10 mm data show highly detailed XAFS well beyond the edge, and temperature variation. The reported in (a) and (c) is the standard error of three repeated measurements, < 0.1% for the 10 mm foil. All XAFS is lost in the 50 mm data, where I 1 /I 0 ' exp(À10).

Dark current characterization
The current amplifier provides the option to null the offset voltage; but doing so makes one blind to the actual value and its variation if negative, and destroys the detector linearity. So, we set the current offset so that the values are always within range; in principle, this does not compromise the sensitivity of the measurement. The resulting nonzero values are referred to as a dark current, the current measured in the absence of X-rays. This was measured throughout the experiment by taking readings at regular intervals with a physical shutter blocking off all incident X-rays (Fig. 4). The counts from the X-ray source are isolated by subtracting the dark-current counts from the total in each of the ion chambers. The intensity ratio due to the change in X-ray flux hence then becomes There are three distinct regions, corresponding to different amplifier tunings and to changing coupled background electronic noise from the beamline. An average or linear model [A + B(t À t 0 )] is fitted in each region. We observe any changes of the dark current by monitoring it through the experiment. The effect of this correction is largest just above the edge, where the dark current makes a larger contribution to downstream ion chamber readings. Here, it increases ln(I u /I d ) by 8.7% for the 10 mm foil, while away from the edge at 11 keV this is around 3%, which continues to be significant. The change in the spectra for the 50 mm foil (Fig. 5) reveals that useful information can be extracted at a high energy, and reveals the power of accurate characterization of the dark current. The standard error of the normalized signal with darkcurrent correction [Fig. 4(d)] is where I represents the counts in each ion chamber and D is the corresponding dark-current count; '0' refers to the upstream detector, '1' refers to the downstream detector. Each is calculated as a standard error. In Fig. 5 the glitch in the uncorrected data at 10.1 keV lines up with the multiple Bragg reflection in Fig. 3(c). This disappears after dark-current correction, as the detectors then have a higher linearity. The numerical corrections applied are summarized in Table 1. It is clear that: (1) dark-current monitoring and correction can be successfully accomplished with a cryostat in low-temperature studies in transmission; and (2) the corrections are very important for thin and thick samples.

Blank measurements and correction for the beam optic
The beam is attenuated by everything in its path -windows into and out of the ion chambers, the cryostat, and air path attenuation. The total attenuation is  Dark-current measurements of each ion chamber (a) I 0 , (b) I 1 , (c) I 2 . The dark current can then be subtracted. The discontinuity follows a change in detector gain. (d) This noise is included in the variance in counts in the uncertainty of the EXAFS spectra. Table 1 Dark current dependence upon time.
The dark-current measurements in Fig. 4 show three distinct regions: 15-30 h; 30-50 h; remainder. I represents the total counts at each ion chamber while sd and se are the standard deviation and standard error, respectively. In the first region of ion chambers I 1 and I 2 , a linear model [D i = A + B(t À t 0 )] is applied.
where subscript s represents attenuation due to the sample and b (blank) represents attenuation by everything else in the beam path between the two ion chambers. We take measurements without any sample to find the absorption of the integrated path. This blank measurement was completed at each energy, so the correction can be applied and experimental errors propagated in a point-wise manner, The uncertainty is then (Chantler et al., 2012a) The 'blank measurements' are energy dependent and depend on the relative gain of the ion chambers. In Fig. 6, at the edge, the systematic effect is in the negative direction, which is opposite to that of the dark current. Removing this effect from the data reduces the attenuation coefficient by up to 52.9% in the pre-edge region, while at higher energies it is increased by up to 3.84%. At the attenuation peak these effects are of similar magnitude, but the combination makes a change to the shape of the distribution. Correction for dark current and blanks contribute to the atomic background of mass attenuation spectra and also to the XAFS oscillations, significantly increasing their sharpness. Changes to the XANES structure and to the k < 3 XAFS are very important (Fig. 7). Hence any theory to understand the XANES, early XAFS or pre-edge should include blank and dark-current monitoring and correction.
The typical accuracy of all measurements to this point are below 0.03%, rising to 0.1% at the absorption edge. It is clear that: (1) dark current and blank monitoring and correction can be successfully accomplished with a cryostat in low-temperature studies in transmission; and (2) the corrections are very important for thin and thick samples.

Additional systematic corrections and the hybrid technique
The monitoring of dark current and blank measurement can lead to large structural corrections to the data (Fig. 9). The measurements of sample, blank and dark current are intrinsic Change in spectra due to blank correction at 10 K for the 10 mm foil. Red: corrected for dark current; black: corrected for blank and dark current, with a reduction of 52.9% pre-edge, and significant change to the slope.

Figure 5
The effect of correcting for dark current in the 10 K data for the (a) 10 mm foil and (b) 50 mm foil. At 10 mm, the largest displacement is near the absorption edge. The largest glitch in the uncorrected data, near to 10.1 keV, disappears after correction for dark current. For 50 mm, correcting for the dark current reveals some consistency in the spectra at high k. It is quite important to monitor the dark current through a spectral data collection.
to both XERT and hybrid techniques. Other systematics that can be very significant are detailed in Appendix A, including: harmonic content of the beam on the sample; energy calibration; thickness and integrated column density determination, and hence the absolute value for the attenuation coefficient; scattering and fluorescence, which yielded a significant correction; roughness of the samples; and bandwidth of the synchrotron beam on the sample. Because the sample was small and could not be mapped, the absolute value of thickness, integrated column density and hence ½= could not be directly measured. However, the nominal or supplier thickness is often quite inadequate. Hence the hybrid technique here uses an independent reference measurement where XERT has been applied and, at least at 298 K, can be directly transferred to normalize and calibrate the absolute value in transmission. Similarly, XERT should measure the energy of the beam at the sample independently and directly. However it is important to realize that a direct measurement of beam energy should be made to enable absolute measurements to be made. But the beamline had no independent measure of energy, so we used, with some limitations of accuracy, selected reference foils with some provenance and the internal edge measurement from the data to calibrate the energy. The measurement of systematic errors requires multiple samples (thicknesses or concentrations) in XERT to calibrate harmonics, fluorescence, roughness and bandwidth; however, this experiment had only two samples and effectively one high-accuracy sample; so estimates and characterization of this, using the hybrid technique, followed a separate XERT measurement on a reference sample. These were all characterized fairly well, with small uncertainty, but indirectly and with some increase in uncertainty as stated. One advantage of the hybrid technique presented here is that it gives a direct insight into systematics and magnitudes for the corresponding, simultaneous fluorescence measurement, discussed later. This is also part of the hybrid technique -to investigate and characterize fluorescence spectra.

Transmission spectra and uncertainty contributions
Fig . 10 provides the X-ray mass attenuation coefficients for the temperature series including roomtemperature data. We were able to investigate or measure several important systematics which can strongly affect many different parameters of the experiment. These include; the energy E or k axis, the amplitude, attenuation or photoabsorption ½= axis, and the structure of the XANES and XAFS   Changes to the 10 K absorption spectra, for the 10 mm foil, due to dark current and blank corrections. (a) Changes to the shape of the slope, pre-edge and edge-jump. (b) The near-edge structure shows significant change in form for processed data with sharper peaks.
oscillations, including dark-current monitoring, blank monitoring, harmonic, thickness calibration, fluorescence, roughness and bandwidth. Importantly, we investigated these for measurements in a cryostat in transmission even though we cannot in this case use a full XERT approach and therefore are using a hybrid approach. There is a large significant difference between the 150 K and 298 K room-temperature spectra, as hoped. All low-temperature spectra show additional resolved peaks at 9.86 keV and 9.94 keV. For 10 K, 100 K and 150 K, the data extend above the zinc K-edge to 11.5 keV; for 50 K, the data extend to 10.1 keV. The error contributions in the absorption spectra are summarized in Table 2. The key contributors to the uncertainty are the counting statistics and the uncertainty in foil thickness. Although the numerous contributions we have discussed produce significant systematic effects on the mass attenuation coefficients, they are well    characterized and so after correction they contribute very little to the total uncertainty. The typical uncertainty in the mass attenuation coefficients is 0.033%.

Background subtraction
The XAFS can be isolated from the background, 0 is the smooth atomic background; Á is the change in absorption from before the edge to the peak of the atomic background, also known as the edge jump; E 0 is the energy of the absorption edge; m e is the mass of the electron; and his Planck's constant.
The software mu2chi (Schalken & Chantler, 2018) follows the common method of fitting a segmented cubic spline; however, it has the advantage of avoiding interpolation of the data onto a uniform k-spaced grid. Hence no information is lost and it also propagates the experimentally derived pointwise uncertainties to allow a properly weighted fit. Fig. 11 demonstrates well resolved sharp peaks in being magnified by over an order of magnitude for the 10 K compared with the 298 K data. The low-temperature data resolves peaks all the way to k ' 17 Å À1 compared with 13.5 Å À1 at room temperature. The broadening at room temperature obscures peaks at k = 7, 9.5 and 11.5 Å À1 .  Table 2 Contributions to the overall uncertainty from each source.

Nanostructure
Corrections are labelled ½= rel if they contribute to the shape of the spectra with an independent contribution to uncertainty at each point. ½= abs notes corrections that relate to broad scale corrections such as the material thickness. The table gives the maximum individual uncertainty and the average contribution . Extra sources of systematic and uncertainty arising from bandwidth and nanoroughness cannot be measured with a single foil in the current setup.  (a) Temperature series of (k). The 50 K data set was collected with a much narrower k-range, up to k = 11 Å À1 . (b) Temperature series of k 2 (k) showing the structural consistency and yet the change in sharpness with temperature, with resolved peaks to k of 17 Å À1 for low temperatures.
we use eFEFFit (Smale et al., 2006;Schalken & Chantler, 2018), which takes the contributions of each potential scattering path theoretically calculated by FEFF and fits them to the XAFS equation, following the popular iFEFFit but with propagated uncertainty. Key parameters are: coordination number (N j ); amplitude reduction factor (S 2 0 ); disorder captured in the Debye-Waller factor ( 2 j ); and effective path distance (r j ), which can be parameterized by a scaling factor j .
The coordination number and inner shell path distances are defined by zinc hexagonal close packing (Ledbetter, 1977). Beyond the nearest dominant coordination shells, a single correction factor can be sufficient. 2 j is estimated from educated guesses about the relations between paths, to reduce the total number of fitting parameters. These follow guidelines by Ravel (2016) and Hudson et al. (1995) and involve trialling different groupings of path parameters in the XAFS equation to minimize 2 r . The transition structure changes at each temperature in the temperature series, especially due to anisotropic thermal expansion and, in other materials, phase changes. Nuss et al. (2010) measured the change in lattice parameters with temperature using X-ray diffraction and found that the expansion was not isotropic; in the hexagonal plane it did not contract below 150 K. Since a is stable at low temperatures (Nuss et al., 2010), we investigated varying the out-of-plane axis c. We implement a grid search through c /a (out-of-plane / in-plane ratio) (Fig. 12). This provides a measure for the lattice spacing and ratio coinciding with a minimum of the 2 r best-fit parameter.
A model based on the 10 K data set analysis was used to provide initial estimates to higher-temperature data sets to minimize the correlation of parameters. A k-range of 4.5-19 Å À1 was used as optimal. Fitting below this range led to a significant increase in 2 r and introduces significant discrepancies with oscillations at higher k -that is, the theory was discrepant below this range. Above 19 Å À1 , or 17 Å À1 for 298 K, the uncertainties dominated over the oscillations. The shortest 37 independent scattering paths (two-leg and threeleg) were modelled. Five scaling and thermal correlation parameters were used in the fit, one for each of the first three independent shortest single leg paths, another to the next 19 longer (r eff < 5.32 Å ) and one for the remaining 15 (5.32 Å < r eff < 7 Å ) ( Fig. 13 and Table 3). The many-body reduction factor should be smaller than 1.0, yet, if left to fit freely, correlations lead to it exceeding 1.15. Fixing S 2 0 = 0.9 following Ekanayake et al. (2021a) was tested, but did not result in any significant change in c /a or fitting parameters.
Each increasing temperature has a smaller magnitude of oscillations and less detailed structure, which reduces 2 r (Fig. 13). At lower temperatures, some fitted peaks are not as resolved as the experiment, and some of the sharper peaks are not resolved which can be due to fitted j being too large. Parameter uncertainties are similar. The increasing sharpness of structures at low temperatures yields deviations of greater significance from the fit, while the structures and fits are qualitatively similar.

Fluorescence: experimental
Fluorescence data were taken using a 100-pixel monolithic germanium detector with a total active area of 5 cm Â 5 cm and depth 7 cm. The beamline software and detector recorded three channels of data: an input count rate (ICR, fast), output count rate (OCR, slow) and integrated region of interest (ROI) (Fig. 14). The third channel only counts photons within the desired energy range around 8.627 keV which is the energy of Zn K. The other two channels are used to correct for detector non-linearities. In fluorescence mode, the relative attenuation coefficient does not correct for systematic effects,

Pixel dead-time corrections
Pixel dead-times are quoted as 0.83 ms and 4.73 ms for input and output count rates, respectively. This XAFS facility does not record the full spectra from the fluorescence detector, in part because of data transfer times and Results of a grid search of c / a in steps of 0.003 with eFEFFit yielding smooth quadratics (Table 3). binary storage space. Input and output count rates (ICR and OCR, respectively) can be used to correct the integrated ROI counts, Non-linearity in pixel response will result in dampening of the XAFS spectrum. Compared with the uncorrected data for a typical pixel, the magnitude is increased by over 20%. Unlike in well characterized absorption data, with matching detectors with high linearity corrected for dark and blank readings etc., the different structure and properties of the fluorescence detector, and limited linearity ranges, means the Bragg peaks (glitches) in intensity do not normalize and cancel out; hence these individual or sometimes two neighbouring data points will have to be deleted during further processing.

Advances in self-absorption correction
There are two major sources of systematic error when performing XAS experiments in fluorescence mode. The spectra have a rising gradient (Fig. 15) instead of a negative gradient which would correspond with absorption coefficients in transmission, or the measurement of absorption. For multipixel detectors, a divergence in the pixels occurs due to the dependence of self-absorption in the sample on the detector geometry. As a photon penetrates the sample, it is absorbed at some depth, and the lower-energy fluorescence photon is Results of eFEFFit modelling at (a) 10 K, (b) 50 K, (c) 100 K, (d) 150 K, (e) 298 K using initial geometry from Nuss et al. (2010). The FEFF fits match the structure very well across a large k-range and also in the lower k-range. At increasing temperature, the oscillations are of smaller magnitude and less detailed. Table 3 eFEFFit fits.
1, 2, 3 and 2 1;2;3 were assigned to the three nearest-neighbour single-scattering paths. 4 and 2 4 were assigned to other paths out to r eff = 5.32 Å , and the rest of the paths out to r eff = 7 Å were assigned 5 and 2 5 . r 1, 2, 3 are the three nearest-neighbour atomic distances defined by the respective . See text.  (2) 0.0058 (2) 0.020 (1) 2 5 (Å 2 ) 0.0054 (6) 0.004 (1) 0.0081 (7) 0.0091 (7) (2) 3.309 (9) r 3 (Å ) 3.8912 (4) 3.912 (7) 3.9079 (4) 3.9153 (4) 4.056 (1) † The 50 K data were collected with a much narrower k-range, up to k = 11 Å À1 . ‡ c = a from Nuss et al. (2010) only included experimental values down to 50 K, hence this is a linearly extrapolated value. § S 2 0 would drift above 1.15 to attempt to accommodate the extra sharp peaks. Hence it was fixed at 1; the fitting compensated by reducing absorbed again by the sample before making it to the detector. Photons that penetrate deeper then have a higher chance of the fluorescence photon being absorbed within the sample, which creates a gradient across the detector. At low attenuation coefficients, e.g. as one goes to higher energies, the photons penetrate deeper into the sample, and a divergence in multipixel spectra is often observed (Trevorah et al., 2019). An example for the ratio of the dead-time-corrected fluorescence signal to the incident count is given in Fig. 15.
Other experimental conditions -such as the use of thin foils, grazing or low angle of incidence, near-90 emission angle, and distance between sample and detector -can help to mitigate these pixel-dependent effects but inevitably the count rate is reduced. To correct for these effects (Chantler et al., 2012b), the program SeAFFluX (Trevorah et al., 2019) was adapted to suit the relevant experimental geometry and the 100-element Ge detector.
The code takes in an experimental geometry and uses FFAST theoretical calculations to fit the self-absorption functional equation (9), Here, inc and out are the angles of the incident X-rays and outgoing fluorescence, respectively. [ f /] is the mass attenuation coefficient of the material at the energy of the fluorescence photon. The penetration depth into the material is more relevant than the material thickness. f is the fluorescence yield, the fraction of events leading to a fluorescent photon emission. Fig. 16 suggests that conditions avoided major divergence of pixel spectra seen in previous fluorescence experiments (Trevorah et al., 2019) -that is, the fluorescence detector was a significant longer distance from the cryostat. The SeAFFluX code, adapted for this experimental setup, has corrected the gradient of the spectra very well. The obvious discrepancy of the two sets is that the magnitude of the XAFS oscillation differs by a factor of $ 3 (Fig. 17). Fluorescence-detector ROIs exclude photons below the desired edge, so that the baseline scaling of the structure can be misaligned and hence mis-scaled. Equation (9) is highly affected by strong oscillations compared with a background or spline scaling, so for example a strong peak (here, ln ratio ' 3.7) is much more heavily damped and a strong valley (here, ln ratio ' 2.7) is (much) less damped. This variation depends upon the emission angle, the log ratio of attenuation and therefore particularly the oscillation magnitude. Here, SeAFFluX uses reference atomic theory, which translates the baseline very well but does not correspondingly scale the peaks and troughs which have strong differential scaling amplification relative to the baseline. So while the base and spline background level is well corrected here (ln ratio ' 3 towards ' 2 above the edge), the data and amplitudes should be used for the amplification of the oscillations -but in selfabsorption these are heavily damped so that the signatures for this in the data are weak. To compensate for the low-magnitude XAFS, a cubic function was fitted through the spectrum and deviations from it are amplified to match the magnitude of transmission XAFS (Fig. 18). As the oscillations are amplified, the noise is correspondingly amplified.

Fluorescence spectra and uncertainties
Figs. 19 and 20 demonstrate results from processing the fluorescence data which produce a very high level of information content, comparable with that in transmission. The typical uncertainty in the fluorescence spectra in ½= Ã is between 0.13% and 0.18%, comparing favourably with the transmission data [see Figs. 11 and 9 from (0.6-1.0) /300 ! 0.25 / 170 ' (0.2-0.3)% and hence ! 0.15% in ½=; and correspondingly 0.0003 -0.001 -0.0015 in ]. The uncertainties are generated by repeated measurements and point-wise propagation by adding uncertainties of dead-time measurements in quadrature. The uncertainties match the level of noise (Fig. 21), confirming that the experimental uncertainties are accurate, to within a factor of two.

Evolution of nanostructure from fluorescence data sets
We use mu2chi to transform the corrected normalized fluorescence spectra in ½= Ã versus E into the XAFS (k). However, the normalization of the oscillations is dependent upon the edge jump Á and the ln ratio in equation (5). Fluorescence measurements are unreliable (zero) pre-edge due to the ROI, so there is no defined edge jump. We include the pre-edge region of the transmission data set to give the appropriate edge jump ratio. Since the fitting procedure only applies above the edge, this provides an accurate amplitude for (Fig. 22). Theoretical fits are then in close agreement with experiment. A strong endorsement for the physical relevance of the data and fit is the continued agreement in the k-range below the Hanning window at 4.5 Å À1 , even down to 2 Å À1 .
The model for each spectrum is identical to that in transmission. 2 r is significantly lower for fluorescence, 4.75 < 2 r < 15.6 (Table 4), compared with Table 3 for transmission, of 11.9 < 2 r < 57. Per se, this is not an argument for neither data set nor fit. Most j and j values are consistent within the uncertainties in both transmission and fluorescence. 1, 4 , 1, 4 correspond to single-scattering paths entirely within the hexagonal plane and 2, 3 , 2, 3 have an out-of-plane component. j and j increase with path length and the j increases with temperature. The asymmetry in j between in-plane and out-of-plane scattering paths is corroborated by Rae et al. (2010). The key difference between transmission and fluorescence is 2 r , which commends both spectroscopies. Near-edge region comparing the transmission spectra with data collected in fluorescence mode and processed by SeAFFluX. The form matches very closely but the magnitude of the XAFS in fluorescence is $ 0.3 that of transmission.

Figure 16
The spread of pixel intensities for the (a) 10 mm foil and (b) 50 mm foil across the whole above-edge energy range. Fluorescence intensities in each of the 100 detector elements, normalized by the upstream ion chamber counts at the absorption edge (9.671 keV), normalized to pixel #55 (close to the centre), across a broad range of energies. For the 10 mm foil, the spread is due to the low count rate and noise, i.e. the experimental uncertainty. For the 50 mm foil, there is a very clear divergence between pixels, caused by self-absorption. The intensity decreases as we move across the detector, as predicted by self-absorption (Fig. 14). SeAFFlux provides geometrical, self-absorption and mode (fluorescence, attenuation) corrections to convert fluorescence intensity into a quantitative ½= Ã pe , ready for analysis using traditional XAFS packages. SeAFFluX uses the experimental geometry to model the divergence and correct accordingly for each pixel, rather than just averaging the pixels. XAS and the XRD results. Both transmission and fluorescence confirm that the single-scattering path 1 is in agreement with XRD results (Nuss et al., 2010), showing little change in lattice spacing in the xy-plane with temperature (Fig. 23).
j has positive trends with temperature as expected from thermal behaviour, except for the 50 K dataset, which is compromised by the short range. Single-scattering paths which have a z-component are consistently larger than in-plane scattering paths, suggesting asymmetric dynamical motion within the crystal. 5 is larger than 4 since, on average, longer scattering paths means larger variance in the variation of the instantaneous path length from thermal motion. The energy offset is slightly large, though it is correlated with other parameters and the Hanning window. 4 is smaller than other scattering paths noting that some of its major contributors are different permutations of scattering between the absorbing atom and nearest neighbours. Theoretical calculations with a cumulant expansion of the XAFS equation (Van Hung et al., 2017) have predicted that the expansion coeffi-cient should quickly go to zero below 100 K, and that 2 should become steady in the range, as observed.
The minima in 2 r follow smooth quadratic curves (Fig. 12), yet show a consistent offset from the literature. The differences in c /a between this work and Nuss et al. (2010) in order of temperature are as follows. For the transmission data sets (T): |Á(c /a) T | = 0.00695, 0.00630, 0.0107, 0.0132 and 0.00724; for the fluorescence data sets (F): |Á(c /a) F | = 0.0140, 0.0143, 0.0117 and 0.0122. In terms of percentages these are jÁðc=aÞ T jð%Þ = 0.382%, 0.345%, 0.588%, 0.720% and 0.391%; jÁðc=aÞ F jð%Þ = 0.770%, 0.787%, 0.643% and 0.665%. XRD measures the distance between two atomic (lattice) sites in the unit cell, whereas XAFS measures the expectation value for the distance between two atoms at any instant. These can be summarized for any single pair of atoms as d 12;XAFS ¼ hr 2 À r 1 i ; where r 1 and r 2 are the vector locations of two atoms and d 12 is the scalar distance as defined either by XRD or XAFS. With any thermal motion for a single-scattering path we therefore expect d 12, XAFS ! d 12, XRD and the difference to be proportional to the magnitude of thermal vibrations. Extended cooling can lead to changes of structure or twinning, and hence a change of the contribution of static disorder, yet none were significant enough to be clearly observed.

Debye model
Comparisons between first-and third-shell Debye-Waller factors ( 2 1 ) are shown in Fig. 24 in which the Debye and Einstein temperatures are fitted to the data. Here we compare with an extended quantum anharmonic correlated Einstein (EQACE) model (Tien, 2021) and with Beni-Platzmann theory (Beni & Platzman, 1976) which treat the vibrational density of states (VDOS) as a delta function at the Einstein frequency, and as a polynomial with a cut-off at the Debye frequency, respectively. We modify the classical anharmonic correlated Einstein (CACE) model which tends to 0 at 0 K, by adding in quadrature a static offset term, from 2 at 10 K. We see good agreement with both the EQACE and Beni-Platzmann theories, which both return plausible results for their respective Einstein and Debye temperatures for the trans- Processing incorporated into SeAFFluX to amplify oscillations and return an accurate set of mass attenuation coefficients. The left hand column (a, c, e) is 10 mm foil data and the right hand (b, d, f ) is 50 mm foil data from the 10 K spectra. Panels (a) and (b) demonstrate the output of SeAFFluX and the shape of the slope. The 10 mm foil transmission spectrum is plotted (black line) in all panels. Panels (c) and (d) show basic spline fits; (e) and ( f ) show results after scaling the divergence relative to this spline. mission data. The fluorescence data, however, demonstrate exceptionally poor agreement for the first shell with all theories returning unreasonably high values of Â E and Â D . Our results do not immediately discriminate between these theories.
A key limitation of some approaches is the absence of any explicit static disorder term in both the EQACE and Beni-Platzmann models. Whilst our modified CACE model contains a static disorder term, the theory is not particularly robust below room temperature. None of these models fully represents the true VDOS of the system; however, the Beni-Platzmann (correlated Debye) model should be preferable. The characterization of dynamic bonding and thermal asymmetry and evolution is clearly observed in the data and modelling, and further work will characterize this more fully.
The precision and apparent accuracy of these results deserves some discussion. Early work (Crozier et al., 1987) claimed to observe no significant evolution of parameters with temperature below 0.02 Å -that is, below 2 pm. Later work (Trö ger et al., 1994) claimed to determine nearest-neighbour bond distances to within 0.015 Å -that is, to within 1.5 pm. Exciting differential or relative XAFS measurements and their potential sensitivity to structure and function (Pettifer et al., 2005) and sequels suggested the possibility of relative measurement of interatomic spacing evolution to 1 fm or 0.00001 Å using novel differential XAFS; which represented a differential Fourier transform in a non-conventional manner. The values presented therein achieved 0.01 Å or 1 pm and the authors commented that XAFS could, in principle, observe structure to 0.001 Å or 100 fm. They showed, however, that the noise floor could reach towards a femtometre scale.
Structural change hypotheses require a detailed assessment of uncertainty or significance as might follow from a 2 r analysis or a Fischer test. In what follows, on the basis of careful measurement and assessment of systematic effects and corrections, we propagate uncertainties and achieve nearest-neighbour shell uncertainties of 0.0002-0.0004 Å or 0.02 pm (20 fm) for transmission with much bigger uncertainties for the limited 50 K dataset; and 0.0002-0.0020 Å or 0.02-0.20 pm for fluorescence. These are consistent with that level. We have measured and reported c /a values from different shell radii, and have assigned estimated uncertainties to these of 0.0015 or 0.030, corresponding to uncertainties of r of shells of order 0.3 pm or so. More importantly, we have compared this with the literature from XRD; part of the observed discrepancy is due to dynamical bond lengths versus site separation, so that our observed Results from modified SeAFFluX after rescaling and averaging over pixels. From top to bottom: accuracy is approximately as stated. The advantage and advance is that data uncertainties were propagated, and that the improvement with standard XAS analysis is of the order of a factor of perhaps 30 compared with the very interesting differential XAFS work.

Supplementary information
Data reference and deposition are becoming very important in the XAS community. Hence it is highly recommended to have data available for other researchers on other beamlines and with other software for analysis. Hence we supply the four datasets for 10 K, 50 K, 100 K and 150 K in transmission in ½= versus E in two formats suggested for cross-platform and reference work -an ifeffit-like .dat file and an IUCr-like .cif file. We supply the four datasets for 10 K, 50 K, 100 K and 150 K in fluorescence in ½= Ã versus E. We supply the corresponding eight versus k data sets for processing. We additionally supply the 298 K versus k data set for processing.
The files labelled Zn-murho-(10-150)K-transmission.cif include the information on the data collection together with seven columns containing the energy, energy error, attenuation coefficient ½=, absolute error in the attenuation coefficient ½= abs , relative error in the attenuation coefficient Results presented as a temperature series for the (a, b) 10 mm and (c, d) 50 mm foils. These share a very strong correlation with our attenuation results.

Figure 21
(a) A sample of the spectra with uncertainties from the 10 K, 50 mm fluorescence spectra. (b) The percentage uncertainty in the fluorescence spectra of 10 K, 50 mm. the relative error in the photoelectric attenuation coefficient ½= PE;rel post-systematic correction from the transmission mode. Files labelled Zn_murho_(10150)K_fluroescence.dat contain three columns detailing the energy, attenuation coefficients and their associated relative error post-systematic and self-absorption correction from the fluorescence mode. The files labelled as Zn_chi_k_ (10150)K_transmission.dat and Zn_chi_k_(10150)K_fluorescence.dat contain the mu2chi output of each of the transmission and fluorescence files, respectively, and contain the k, , Á and E 0 values that can be fed directly into efeffit for analysis. Finally the file Zn_chi_k_298K_transmission.dat also contains the post-mu2chi data from Ekanayake et al. (2021a,b) which we compare directly with in this work; the original attenuation data files can be found in the supplementary information section of those papers. That is, a total of 17 data sets are supplied in the supporting information. Headers explain columns and format.

Conclusion
We report the first transition metal XAFS using the hybrid technique and at low temperatures. This is also the first (hybrid-like) experiment at the Australian Synchrotron. We Fluorescence XAFS spectra in k 2 versus k and theoretical fits at (a) 10 K, (b) 50 K, (c) 100 K and (d) 150 K. The fit follows the data very closely, even below the fitting window.

Figure 23
Optimized 2 r for c / a from the grid search and compared with the literature (Nuss et al., 2010). An uncertainty of 0.003 follows the error analysis, step size and the uncertainty of the quadratic (Fig. 12, Table 3). Table 4 Fitted parameters using fluorescent data sets; defined as for the transmission measurements,  (5) report high-accuracy measurement of the XAS of zinc metal at the Zn K-edge. We report a methodology for cryogenic and temperature series measurements which is robust, valid for both transmission and fluorescence XAS, and attains parameter accuracy. We demonstrate the use of SeAFFluX as a precision tool for converting fluorescent spectra to extract the absorption XAFS on an absolute scale. We observed thermal evolution and dynamic lattice parameter ratios in very good agreement with literature XRD values, noting the expected difference between dynamic bonding (XAFS) and site separation (XRD). We are able to investigate aspects of thermal behaviour models using this approach. More work is called for.

APPENDIX A Harmonic contributions
The equipment was mounted at a beamline which did not have a pre-mirror between the source and the double-crystal monochromator. Such mirrors are frequently used as a high cut filter to eliminate high-energy radiation which might interfere with the experiment. A Si(111) double-crystal monochromator was used in this experiment. The Si(222) reflection is forbidden; so the major contributor to harmonics is the Si(333) reflection at an energy three times that of the fundamental energy.
Since the third-order harmonics in general have a much lower attenuation coefficient, the linear relationship between thickness t and absorption ð½=½tÞ is not valid. Thicker foils will absorb most of the fundamental energy photons and higher-energy photons will then dominate in the downstream beam. Quantifying this non-linearity gives where is the proportion of harmonic content, the effective harmonic content (Tran et al., 2003;de Jonge et al., 2005;Glover & Chantler, 2009), [/] f is the mass attenuation coefficient at the fundamental energy, and [/] h the mass attenuation coefficient at the harmonic energy. An extra parameter DC offset is added to check the robustness and accuracy of the earlier determination of dark-current measurements -it is zero within uncertainty. Normalization of the local daisy wheel foil thicknesses in the beam is performed at the highest energy, where harmonic contamination is considered negligible; in this case 11.29 keV (Tran et al., 2003). Nine aluminium foils in the daisy wheels were used. Two independent methods may be used to determine the daisy wheel foil thicknesses, relying on either the nominal thicknesses or FFAST attenuation tables: (1) The first uses the experimentally derived lnðI 0 =I 1 Þ = ½=½t and mass attenuation coefficients [/] from FFAST (Chantler, 1995(Chantler, , 2000. The ratio between these two, divided by the known density, retrieves the thickness. (2) The second plots measure lnðI 0 =I 1 Þ = ½=½t against [t] from nominal thickness values and then fits a straight line to determine [/]. Then, t is determined with fitted values of [/]. The two methods differed by < 10 À5 %. Harmonic measurements were taken at various instances during the experiment, spanning the whole energy range. The results presented limit the dark-current offset to about one standard deviation as in Table 1.
With most energy points returning indistinguishable from 0 and others below = 5 Â 10 À6 , and given that the intensity downstream with no foils is between 2 Â 10 5 and 3 Â 10 5 counts, this equates to a maximum contamination of 1-2 counts at each point. The dark-current standard deviation is $ 1.5-2.5 counts in the upstream and downstream ion chambers. By substituting a value for the effective harmonic parameter of 5 Â 10 À6 into equation (11), we observe that any harmonic contamination will be negligible. Hence we were able to investigate this carefully, but it was not a major source of systematic error in this experiment.

APPENDIX B Energy calibration
Hysteresis in the positioning motors of the monochromator can cause the energy to drift over large energy ranges; hence an independent energy measurement is desired. Standard XAFS uses a single calibrated reference foil. Used in its original form, XERT should directly measure the energy of the beam through the sample with an independent powder or crystal reference. In the absence of this at this beamline, we used zinc, copper and tantalum reference foils (Fig. 25). The edge energy in a material is most commonly defined experimentally as the first peak in the derivative spectra (Kraft et al., 1996). Energy resolution and photon bandwidth affect this energy in an ill-defined manner. Using this definition, the Zn K-edge was 9663.23 (18) eV, the Cu K-edge was 8981.87 (21) eV and the Ta L 3 -edge was 9880.13 (10) eV, compared with reference values of 9660.755 (30) eV, 8980.476 (20) eV (Kraft et al., 1996) and 9881.1 eV (Wong, 1999). 1 The zinc K-edge energy is also determined at the four temperatures from the collected data sets, which have a much finer point spacing than the reference foil sets.
In Fig. 26 a Gaussian has been fitted with a linear background to deal with asymmetry. Derivative peak values ranged from 9661.23 eV to 9.661.30 eV with standard deviations < 0.1 eV, which is a consistent offset from the value 9660.755 (30) eV (Kraft et al., 1996) of 0.5 eV. The consistency for all temperatures suggests that the monochromator was consistent over time, and may correspond to a single energy offset of 0.5 eV (Table 5). This correction is important and useful, but not that definitive and subject to further investi-  gation in the future, especially given the mixed evidence provided by Table 6. However, one can argue that the XAFS is simply offset by this amount and that XAFS fitting software is therefore somewhat independent of this correction.

APPENDIX C Thickness determination
To convert from ½=½t to ½= (cm 2 g À1 ) requires a precise measurement of the foil thickness at the point where it interacts with the X-ray beam; or scaling to a calibrated reference. A recent experiment on room-temperature zinc foils made an accurate measurement of foil thickness allowing scaling of these ½=½t spectra (Sier et al., 2022; Ekanayake et al., 2021a,b). The column density was determined by selecting a small well defined energy range in the pre-edge spectra of the 150 K and room-temperature 298 K data, 1 ½t 150 K ðcm 2 g À1 Þ ¼ ½= 298 K ðcm 2 g À1 Þ = ½ ½t 150 K : The scaling factor calculated was 87.372 cm 2 g À1 for a density of 7.112 g cm À3 . The foil is at 45 so corresponds to a foil thickness of 11.380 mm which is within the manufacturer specifications of 10 mm AE 15%. The accuracy in thickness claimed in the room-temperature series is < 0.024%. Measurements of thermal expansion with temperature using XRD have some uncertain error bars (Nuss et al., 2010). The c / a ratio is the dominant variation; from 10 K to 150 K this is 5.5 Â 10 À5 K À1 ; from 150 K to 298 K it is roughly 6.2 Â 10 À5 K À1 , fairly uniform and monotonic. Theoretical calculations suggest that the thermal expansion via a cumulant approach should go rapidly to zero below 100 K; similarly, experimental XRD results ( Determination of the Zn K absorption edge by finding the first peak in the derivative spectra for (a) 10 K, (b) 50 K, (c) 100 K and (d) 150 K. The energy is determined to a fraction of the point spacing.

Table 5
Determination of the absorption edge by finding the first peak in the derivative spectra.
Each temperature is self-consistent to within error and demonstrates a potential calibration error of 0.5 eV.  † Kraft et al. (1996). Table 6 Results from fitting the reference foils are inconsistent -reference and deviations do not follow a physical progression with energy. The nanoroughness of the sample is a change in the integrated column density [t] through the beam footprint, including surface roughness and internal roughness. It affects the accuracy, especially for thin foils like those we are using (de Jonge et al., 2007). The change in the mass attenuation coefficient due to nanoroughness is  Effect of correcting for the scattering according to equation (13). The model is effective across the whole energy range. 'XERT' is the reference data taken at 298 K (Ekanayake et al., 2021a,b;Sier et al., 2022). '150 K' is this work, for T = 150K.
measurements conducted over a range of different sample thicknesses, we cannot independently investigate the sample nanoroughness with these data. Given that this experiment used the same sample foil as Ekanayake et al. (2021b), we take the value obtained in that analysis, [t] = (9.9 Â 10 À4 AE 4 Â 10 À6 ) g cm À2 for the 10 mm sample, with a 10% uncertainty. This results in a correction of the data of up to 2.52% for the 10 mm sample. The effect of the energy bandwidth in the incoming X-rays, on the spectra, should also be considered. The Si(111) monochromator allows for a small bandwidth of energies that are centred around our desired E 0 due to intrinsic bandwidth, heat load and thermal and mechanical stresses within the monochromator. The absorption can be modelled by where I represents the intensity around an infinitesimal energy range dE centred on E 0 and ½= m is the measured mass attenuation coefficient at nominal energy. The effect of this will be seen predominantly as a broadening in areas of high absorption gradient such as an absorption edge, and should also grow with sample thickness.
With only a single foil thickness, we can take a similar approach as with the nanoroughness, and compare with earlier calibrated data (Ekanayake et al., 2021a,b) to observe the signature. In comparing two temperatures (150 K and room temperature), we expect to see some thermal broadening. By looking at the sharp first edge of the peak, we see where the thermal broadening is at its lowest, and the bandwidth effect is at its highest. Still in this range, the thermal broadening is dominant with no evidence of further effects from bandwidth, hence we can progress confidently without significant bandwidth effects.